Correlation energies in the random phase approximation using realistic interactions 
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The self-consistent random phase approximation (RPA) based on a correlated realistic nucleon-nucleon in- 
teraction is used to evaluate correlation energies in closed-shell nuclei beyond the Hartree-Fock level. The 
relevance of contributions associated with charge exchange excitations as well as the necessity to correct for the 
double counting of the second order contribution to the RPA ring summation are emphasized. Once these effects 
are properly accounted for, the RPA ring summation provides an efficient tool to assess the impact of long-range 
correlations on binding energies throughout the whole nuclear chart, which is of particular importance when 
starting from realistic interactions. 
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On the mean-field level, the nuclear many-body problem 
can be treated in the well known Hartree-Fock (HF) ap- 
proximation. In conjunction with phenomenological inter- 
actions, this scheme is successful in describing various nu- 
clear ground state properties. Recently, realistic nucleon- 
nucleon (NN) interactions have also been regulated for imple- 
mentation in nuclear structure calculations in two novel ap- 
proaches, (i) low-momentum NN interaction from the renor- 
malization group theory, Viow* JjJ] , and (ii) correlated interac- 
tions Vucom constructed in the framework of unitary correla- 
tion operator method |2|. 

In connection with these interactions, the HF approxima- 
tion underestimates the binding energy due to the inadequate 
representation of long-range correlations. In recent HF stud- 
ies based on Vi w* and Vucom, these missing correlations 
have been recovered within many-body perturbation theory 
(MBPT) by evaluating corrections up to third order |3, |4[ . In 
general, correlations beyond HF are more relevant in studies 
based on realistic NN interactions than in phenomenological 
models based on, e.g., Skyrme or Gogny functionals, because 
the latter already mimic part of the many-body correlations 
through the phenomenological fit. Conversely, studies based 
on realistic interactions allow for more profound insights into 
the many-body dynamics 0. 

Going beyond low-order MBPT, correlation energies can 
be obtained using the random phase approximation (RPA) as 
a tool to evaluate a partial summation over particle-hole ring 
diagrams. In the language of RPA, ground state correlations 
emerge from a coupling to giant resonances and surface vibra- 
tions | g, |7|, |8|, |9[ • The correlation effects of these two classes 
of collective motion are different — giant resonances influence 
the binding energies, while surface modes have a more pro- 
nounced effect on charge density distributions flfll fTHl . By 
applying the RPA method, these excitations are approximated 
as harmonic vibrations. The total binding energy can then be 
evaluated as the sum of the zero-point energies of all the possi- 
ble modes. The RPA framework has been employed in several 
studies of the ground state correlation energies mainly associ- 



ated with quadrupole and octupole modes, and pairing vibra- 
tions Ofmfl^TTlllrill . It should be noted that for systems 
made of more than one fermion species, other modes, e.g. 
Gamow-Teller transitions in nuclei, become possible. How- 
ever, no discussion is usually found in the literature regarding 
the relevance of charge exchange excitations. 

In this work we study RPA correlation energies in con- 
junction with correlated realistic NN interactions derived from 
the Argonne V18 potential 1 17] in the framework of the uni- 
tary correlation operator method (UCOM). The short-range 
central and tensor correlations induced by the realistic poten- 
tial are described by a unitary state-independent transforma- 
tion 0, fl^l . The central correlations induced by the strong 
short-range repulsion (core) of the interaction are described 
by a unitary shift in the relative coordinate of all nucleon 
pairs Q. [T3l . The tensor part of the interaction induces strong 
correlations between spin and relative orientation of the nu- 
cleons. Only the short-range system-independent part of these 
tensor correlations is described explicitly by the unitary trans- 
formation, the effect of long-range tensor correlations as well 
as all residual correlations has to be described by the many- 
body method employed. 

In this framework we obtain a correlated interaction Vucom 
which is phase-shift equivalent to, but much softer than the 
original potential. Technically, the unitary transformation 
modifies the off-shell behavior of the interaction by introduc- 
ing strong non-local terms. The application of Vucom in no- 
core shell model calculations for 3 H and 4 He shows a dramatic 
improvement of the convergence behavior flill . At the same 
time, the tensor correlator can be tuned in order to minimize 
the effect of the net three-body force, which is the sum of the 
genuine three-body potential accompanying the bare NN po- 
tential and the three-body terms induced by the unitary trans- 
formation of the Hamiltonian. As no-core shell model calcula- 
tions show, the two-body Vucom provides a good quantitative 
description of ground and low-lying excited states throughout 
the p-shell ll9ll . Using the same correlated interaction Vucom 
we have performed HF and RPA calculations for closed shell 
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nuclei throughout the nuclear chart Based on the HF 

solution the effect of long-range correlations on the binding 
energies was estimated via MBPT and good agreement with 
experimental binding energies was found |4]. 

In the present work we investigate the effect of long- 
range correlations in the framework of a fully self-consistent 
RPA 1 20]. Using the single-particle basis resulting from HF, 
the RPA configuration space is built and the generalized eigen- 
value problem posed by the RPA equations is solved, 
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where the eigenvalues a> v correspond to RPA excitation ener- 
gies. The forward and backward-going particle-hole ampli- 
tudes, X v ph and Y",, respectively, are related to the transition 
amplitudes to the excited states \*¥ v ) by 



OPvlcfopPo) = S a , p Spj, X v ph + 5 a ,hi 



(2) 



Here and in the following, the indices p (h) always label par- 
ticle (hole) states while Greek letters refer to any orbits. c T (c) 
are the usual creation (annihilation) operators. The HF+RPA 
scheme is applied in a fully self-consistent way, i.e. the same 
translational invariant Hamiltonian H mt = T — T cm + Vucom, 
which formally is a two-body operator, is used in the HF equa- 
tions that determine the single-particle basis and in the RPA 
matrices A and B. This ensures that the RPA amplitudes do 
not contain spurious components associated with the center- 
of-mass translational motion 1 20] . 

There are several approaches to calculate the RPA ground- 
state correlation energy. In the present study, we consider two 
different formulations, based on (A) the direct evaluation of 
the expectation value of the Hamiltonian and (B) the quasi- 
boson approximation. The expectation value of the Hamilto- 
nian Hint (or, in general, the two-body part of the Hamiltonian) 
is given by 
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where H a ^ 7 g are the antisymmetrized matrix elements of H[ nt . 
Eq. is still exact. By approximating the transition am- 
plitudes according to Eq. and employing the RPA equa- 
tions Q one is led to express the total energy of the system as 



E — Ehv + £r 



where the HF binding energy supplemented with the RPA cor- 
relation energy, 
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In the quasi-boson approximation the ground state energy 
is described as the zero point energy of a collection of har- 
monic vibrations. Using the oscillator-projection method by 
Rowe |7] one obtains, 



(6) 
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Eqs. (0 and (|6} employ the RPA eigenvalues a> v and 
backward-going amplitudes Y v h as well as the HF single- 
particle levels e a . From the completeness insertion in the last 
row of Eq. Q and the summation over all (proton and neu- 
tron) orbits, it is clear that one needs to sum over all possible 
final states predicted by the RPA formalism, independently of 
isospin. Hence, the sum over v in Eqs. (0 and (|6jl includes 
not only all the multipolarities and parities but also charge- 
exchange processes. 

In comparison to a order-by-order summation of particle- 
hole ring diagrams in MBPT, the above formulae implicitly 
double count the second order contribution, as pointed out by 
Fukuda et al. |6] and Ellis 12 ill . To show this, we expand 
the RPA eigenvalues a> Y and the amplitudes X ph and Y v ph in a 
perturbation series of the interaction H mi , as done in Ref. |6]. 
Inserting these into Eqs. (jSJi and © one obtains, 
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where £ <2) corresponds to the second order contribution in 
MBPT 
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and £V to the contribution of the ring diagram at third order 
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The double counting of the second order contribution £ <2) is 
evident and has to be corrected for explicitly. Beyond the 
second order, EjJ A does not introduce any further overcount- 



ings while E^ A remains somewhat troublesome. The 
double counting of £ (2) is intrinsic to the quasi-boson approx- 
imation and can be avoided only in a formalism that (beyond 
the RPA approach) explicitly recouples particle and hole states 
between different phonons. This can be achieved at the level 
(4) of the many -body self-energy 
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1: (Color online) Partial contributions to the RPA correlation 
energies in 40 Ca, 90 Zr, and 208 Pb as a function of multipolarity J. 
The correlation energies are evaluated using Eq. |6) (left panels) and 
Eq. J5J (right panels). Solid and dashed lines correspond to contri- 
butions from natural and unnatural parity excitations, respectively. 
Contributions from charge-exchange excitations (C.E.) for natural 
and unnatural parities are represented by dot-dashed and dotted lines. 
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FIG. 2: (Color online) Binding energies per nucleon for a series 
of closed-shell nuclei in comparison with experiment. Shown are 
the HF energy, the energy including RPA correlations according to 
Eqs. O and © (labeled HF+RPA (A) and HF+RPA (B) , resp.), and 
those including RPA correlations corrected for the double counting 
of the second order (denoted HF+RPA (A B) -MBPT). All calculations 
are based on the correlated Argonne VI 8 interaction. 



One has to keep in mind that the particle-hole ring sum- 
mation does not include all possible diagrams, e.g., the third 
order ring term £ (3) is only one of three third order contri- 

° ring J 

butions and two-particle two-hole diagrams are neglected by 
both Eq. and 0. The latter are known to be approxi- 
mately compensated by Pauli exchange effects between dif- 
ferent phonons at high order |'23| • 

For a microscopic theory based on realistic nucleon- 
nucleon interactions, the correlation energy beyond the sim- 
ple HF approximation is sizable. This remains true even 
after regularizing the interaction to account for the effects 
of short-range correlations. The residual long-range corre- 
lations manifest themselves in a sizable second-order contri- 
bution, £ (2) |4J. Therefore, the naive application of Eqs. (|5} 
and l|6) would lead to a strong overestimation of the corre- 
lation energy. In the context of phenomenological models 
Eq. (HJ is often applied without corrections for double count- 
ings 1121 1 1 3l 1 1 4l Il5l llol . This might be compensated by 
also neglecting the contribution from charge exchange terms. 
However, these are two (not well controlled) errors that do not 
necessarily cancel each other 1 24] . 

We apply this scheme to evaluate correlation energies for 
closed-shell nuclei throughout the nuclear chart based on a re- 
alistic NN interaction. We employ the Vtjcom derived from 
the Argonne V18 interaction using the optimal correlation op- 
erators determined in Ref . I18II . The range of the tensor cor- 
relator in the triplet-even channel was fixed to reproduce the 
binding energies of A < 4 nuclei in no-core shell model cal- 
culations (/„ = 0.09 fm 3 ). The same correlator was used suc- 
cessfully in HF and MBPT calculations reported in Ref. |4B. 
For a systematic calculation of the RPA correlations, Eq. {0 
was projected onto good angular momentum 7 and parity n 
and all available multipolarities were taken into account. All 



calculations were performed using 1 3 major harmonic oscilla- 
tor shells. The spurious contributions in the 1" channel were 
excluded from the calculations of correlation energies. 

In Fig.^we display the individual contributions to the cor- 
relation energies evaluated in the two different formulations 
(Eqs. d3 and ©) for 48 Ca, 90 Zr, and 208 Pb. Shown are the 
contributions to E^ A and E^ A as function of the multipolar- 
ity 7* = 0* - 13* separated into natural parity, n = (-l) y , and 
unnatural parity, n = (— excitations as well as charge 
exchange excitations. In general, the RPA correlation energy 
increases with J, reaches the maximum for 7 = 3-4, and 
slowly decreases towards higher multipolarities. Both, natural 
and unnatural parity states are equally important and charge- 
exchange excitations also have significant contributions to the 
correlation energy. Although in all nuclei the largest contribu- 
tions come from the 7 = 3-4 excited states, one obviously 
needs to include all other multipolarities as well. This is es- 
pecially important for heavier nuclei, where the correlation 
energies are widely distributed over various multipolarities up 
to 7= 13. 

The overall sum of correlation energies displayed in Fig.^ 
provides the correction to the binding energies in finite nu- 
clei. In Fig.|2we show binding energies per nucleon for sev- 
eral closed-shell nuclei obtained in HF with and without the 
inclusion of the RPA correlation energies, in comparison to 
the experimental binding energies 1 25 ] . The plain HF calcu- 
lations underestimate the binding energies due to the inade- 
quate description of long-range correlations. Inclusion of the 
correlation energies resulting from Eqs. (|5jl and l|6) without 
correction for the double-counting of the second order con- 
tribution leads to a strong overbinding. As discussed by da 
Providencia [26], only after explicit correction for the dou- 
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the correlation energy. Then the RPA ring summation pro- 
vides an efficient tool to evaluate correlation energies through- 
out the nuclear chart. In connection with Vucom, the RPA cor- 
relation energies generally confirm the results of second or- 
der MBPT, indicating that the net contribution of higher-order 
ring-diagrams are moderate although the second order itself 
is large. Both, RPA ring summation and low-order MBPT, 
provide efficient tools for nuclear structure calculations with 
correlated realistic NN-interactions. 

This work is supported by the Deutsche Forschungsgemein- 
schaft (DFG) under contract SFB 634. 



FIG. 3: (Color online) Binding energies per nucleon including 
RPA correlations and second-order correction (HF+RPA <A B) -MBPT) 
compared to energies resulting from second-order perturbation the- 
ory only (labeled HF+MBPT). Compare to Fig.|2] 

ble counting by subtracting the second order contribution, i.e. 
by using £hf + £rpa - £ <2 \ we obtain a proper estimate for 
the binding energy including correlation effects which is in 
good agreement with experiment. The difference between the 
two schemes for the evaluation of the RPA correlation energy 
traces back to the overcountings of higher-order terms in Eq. 
0. Only E ( ^ A resulting from Eq. ijfji corresponds directly to 
the particle-hole ring summation. 

For a more detailed discussion, in Fig. [3] we show the cor- 
rected RPA energies, E HF + £rpa - E <2 \ in comparison to the 
direct second order perturbative estimate, £hf + £ <2) - The 
binding energies per nucleon obtained by employing the sec- 
ond order perturbation theory agree rather well with the re- 
sults of the RPA ring summation. Although the second-order 
correction to the HF energy is large, the higher order contri- 
butions included in the ring summation seem to have a rela- 
tively small net effect. Beyond 48 Ca the binding energies per 
nucleon including E^ A are in excellent agreement with the 
perturbative second-order results and also with experimental 
values. The ring summation provides systematically larger 
correlation energies than the plain second-order. This is in 
line with no-core shell model calculations for light isotopes 
which also predict energies somewhat lower than the second- 
order estimate 1 19]. 

In conclusion, we have employed a self-consistent RPA ap- 
proach to evaluate correlation energies based on a correlated 
realistic NN-potential. Correlations beyond the HF level have 
a sizable impact even if one uses regularized interactions like 
Vucom- We point out the need to sum over all possible ex- 
citation modes, including charge-exchange excitations, and to 
correct for the double counting of the second-order contribu- 
tion, when using standard expressions like Eq. (|6jl to evaluate 
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